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ABSTRACT 

The coefficients defining the mean electromotive force in a Galloway-Proctor flow are deter- 
mined. This flow shows a two-dimensional pattern and is helical. The pattern wobbles in its 
plane. Apart from one exception a circular motion of the flow pattern is assumed. This cor- 
responds to one of the cases considered recently by Courvoisier, Hughes and Tobias (2006, 
Phys. Rev. Lett., 96, 034503). An analytic theory of the a effect and related effects in this flow 
is developed within the second-order correlation approximation and a corresponding fourth- 
order approximation. In the validity range of these approximations there is an a effect but no 
7 effect, or pumping effect. Numerical results obtained with the test-field method, which are 
independent of these approximations, confirm the results for a and show that 7 is in general 
nonzero. Both a and 7 show a complex dependency on the magnetic Reynolds number and 
other parameters that define the flow, that is, amplitude and frequency of the wobbling motion. 
Some results for the magnetic diffusivity rjt and a related quantity are given, too. Finally a re- 
sult for a in the case of a randomly varying flow without the aforementioned circular motion 
is presented. This flow may be a more appropriate model for studying the a effect and related 
effects in flows that are statistical isotropic in a plane. 
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1 INTRODUCTION 

In the astrophysical context, turbulent flows, e.g. in stellar con- 
vection zones or in accretion discs and galaxies, are generally 
anisotropic and time-depe ndent. A simple model of a flow with 
such properties is that by iGallowav & ProctoJ ( Il992h . This flow 
is two-dimensional, depends only on two Cartesian coordinates, 
e.g. X and y, which can simplify the analysis significantly, even 
in dynamo problems that are inherently three-dimensional. The 
Galloway-Proctor (GP) flow is related to a flow considered by 
Roberts (1972). The Roberts flow0 is an early example of a spa- 
tially periodic flow that produces an alpha effect. The alpha term 
in the averaged form of the induction equation is crucial to model 
the generation of large-scale magnetic fields from small-sc ale heli- 
cal flu i d motions in s tars a nd galaxies; see, for exa mple, iMoffatj 
( Il978h . IParkel jl97 9'"). and 'K rause & Ra dlef (1980) for standard 
references. However, unlike the Roberts flow, the GP flow is time- 
dependent with a flow pattern wobbling in the {x,y) plane in a 
circular fashion. Both the GP flow and the Roberts flow have a ve- 
locity component out of this plane such that the flow can be fully 
helical, i.e. the velocity is proportional to its curl. 

Particularly important is the dependence of the a effect on 
the magnetic Reynolds number, Rm- While for the Roberts flow a 
declines with Rm in the large Rm limit, in the case of the GP flow 
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according to the results by ICourvoisier. Hughes & To bias" ( 200^ ) 
(in the following referred to as CHT06) and ICourvoisier, t2008l) 
there is a more complicated dependence on Rm with sign changes 
and no indication of convergence with increasing R^. 

In many studies turbulent astrophysical flows have been mod- 
elled by random forcing. In the case of of helical isotropic turbu- 
lence such investigations show that a approaches a finite value as 
soon as i?m exceeds a value of the order of unit y. This has bee n 
observed at least for Reynolds numbers up to 200 dSur et alj2008h . 

The purpose of this paper is to study the effects of the GP 
flow in more detail in order to understand the influence of time- 
dependence and anisotropy on the value of a and other turbulent 
transport coefficients. In particular, it is important to document the 
differences and similarities with turbulent flows that are statistically 
isotropic and irregular in space and time. We focus attention here 
on the simplest case considered in CHT06 with a flow being purely 
periodic in time and add a few results for a simple flow with random 
time dependence. 

A number of similarities, but also some striking differences 
between turbulent flows and the Roberts flow are known. Similar 
in both flows is the fact that there is an a effect whose magnitude 
increases with _Rm as long as the latter does not exceed some value 
in the order of unity. However, for larger Tim, t he a coefficien t 
in isotropic turbulence settles to a constant value dSur et al]|2008h. 
while for the Roberts flow a ten ds to zero as Rm — > 00 jSowarj 
1 1981 1 19891 : iRiidler et al.l l2002a''b'). Furthermore, there is no 7 ef- 
fect, or pumping effect, neither for isotropic turbulence nor for the 
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Roberts flow. On the other hand, in the time-dependent GP flow 7 
effects have been reported (CHT06). A 7 effect corresponds to an- 
tisymmetric contributions of the a tensor. This raises the question 
about the possible existence of antisymmetric contributions to the 
turbulent magnetic diffusivity tensor, or rjt tensor. 

These aspects are now straightforward to address using the 
recently developed test-field method to calculate numerically all 
components of the ct and rjt tenso rs defining the mean electro mo- 
tive force f for a given flow field jSchrinner et"aDl200l l2007h . If, 
as we assume here, too, the mean magnetic field B depends only 
on one of the Cartesian coordinates, say z, only two 2x2 ten- 
sors for a and ijt are of interest. The test-field method has recently 
bee n used to calculate diag onal and off-diagonal components of 
rjt terandenburg et alj|2008l). the magn etic Reynolds number de- 
pendence of a and rjt jS ur et al."2008'), as well as their scale de- 
pendence terandenburg, Riidler & Schrinner 2008). We begin by 
exploring general properties of the mean electromotive force in 
the GP flow and present analytical results for coefficients like a 
and 7, which are crucial for the electromotive force, gained in 
the second-order correlation approximation and in a corresponding 
fourth-order approximation. After explaining the test-field method 
we give a series of numerical results for such coefficients, which 
are independent of approximations of that kind, and discuss them 
in detail. 



2 MEAN-FIELD ELECTRODYNAMICS WITH 
GALLOWAY-PROCTOR FLOW 

2.1 Definition of the problem 

Consider a magnetic field B in an infinitely extended homogeneous 
conducting fluid with constant magnetic diffusivity 77 moving with 
a velocity u. Its behavior is governed by 

dtB -TjV'^B -V X {ux B) =0, V-B = 0. (1) 

Referring to a Cartesian coordinate system [x, y, z) the velocity u 
is specified by 

u = — i; X V^p — z knip (2) 
with 

= [cos{kiix + tpx) + COs(fcHj/ + 'Py)\ ■ (3) 

Here z means the unit vector in the z direction, fcn is a positive 
constant such that 2n / kn is the length of the diagonal of a flow cell, 
and and ipy are functions of time to be specified later. Further 
we have uq = Urms/V^- In the special case (fx = (Py = the 
flow agrees with a Roberts flow. For non-zero ipx or (fiy a properly 
moving frame of reference can be found in which we have again 
a steady Roberts flow pattern. In our original frame each point of 
this pattern moves with the velocity —k^^{dtfx, dtfy) in the xy 
plane. In (O the ratio of the flow components in the xy plane and 
in z direction has been fixed such that the modulus of the average 
of the kinetic helicity u • (V x u) over all x and y for given uo 
takes its maximum. With the signs chosen this average is equal to 
-2u|fcH. 

In view of the first example treated in CHT06 we specify the 
flow generally defined by ^ and ([3} further to be a Galloway- 
Proctor flow and put 

(fix = t COS ujt , tpy — t sin uit , (4) 



where e and uj are considered as non-negative constants. We label 
this flow in what follows by (i). Each point of this pattern moves 
with the frequency uj/2tt ona circle with the radius e/ku- 

To come closer to a turbulent situation CHT06 added a ran- 
dom function of time to the arguments u)t in Another case of 
some interest occurs if we simply interpret ip^, and ipy as random 
functions. More precisely we put 

= f^<t>x{t/Tc) , ify = €(t>y{t/Tc) , (5) 

where e is again a constant, <j)x and <j)y are two independent but 
statistically equivalent random functions, which take positive and 
negative values between —1 and 1 and tend to zero with growing 
moduli of the argument, and Tc is some correlation time. We label 
this random flow by (ii). 

2.2 Mean-field concept 

Adopting the mean-field concept, we denote mean fields by an 
overbar and define them as averages over all x and y. We have 
then u = 0. Taking the average of (T} we find 

dt'B - rjV'^B - Vx£ = 0, V-B = 0, (6) 

with the mean electromotive force 

£=TI^, (7) 

where b = B — B. [In reduces simply to (0, 0, 9z).] From 

l[T) and l|6) we conclude that b has to obey 

(9t-77V^)6 = (B ■ V)u - (u • V)B 

+V X (u X 6 - uxb) , V ■ 6 = . (8) 

We adopt here the assumption that the mean electromotive 
force £ is, apart from u and rj, completely determined by B and its 
first spatial derivatives. (This assumption will be relaxed in Sect.|4]) 
This implies that there is no small-scale dynamo and that sufficient 
time has elapsed since the initial instant so that £ no longer depends 
on any initial conditions. Since B is by definition independent of 
X and y its spatial derivatives can be represented by V x B. We 
write simply J instead of V x B, being aware that the mean elec- 
tric current density is really V x B/ji (rather than J), where ji 
is magnetic permeability of the conducting fluid. Clearly we have 
now J — {—dBy/dz, dBx/dz, 0). For the sake of simplicity we 
further assume that B is steady. In the so defined framework we 
may write 

£i = ciijBj - rjijjj (9) 

with tensors aij and rjij determined by u and rj only. Both aij and 
rjij, and so £i, too, depend in general on time. 

We see from ([8} that, if B is a uniform field, b and therefore 
£ are independent of Bz- Hence we have = 0. Furthermore, 
since J z ~ 0, clearly rjis is without interest, and we put rjis = 0. 

2.3 Mean electromotive force in case (i) 

For a more detailed investigation of £ we focus on the fluid flow of 
type (i). In this case atj and rjij are periodic in time with a basic 
period equal to that of u, that is 2tt/ui, or (as we will see below) a 
fraction of it. 

Remarkably the velocity field u = u{x,y,t) defined by 
([3) and l|4j is invariant under a 90°rotation about the z axis and a 
simultaneous retarding by n/2ij (that is, cot ^ cut — tv/2). Conse- 
quently the aij in the correspondingly rotated coordinate system, 
which we denote by a'^j, have to satisfy the relation 
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ot'ij {t — 7r/2ti;) = atj (t) . 



(10) 



If we consider for a moment the change of the spatial coordinate 
system only and ignore any time dependence we have a'n — Q22, 
a'i2 = —021, a'21 = —ai2, a'22 = ail, asi = — 032 and 032 = 
031. Hence l IlOl ) provides us with 

aii(i) = +Q22(t - 7r/2aj), a22(t) = +aii(t - 7r/2i^), 
ai2{t) = -a2i{t-Tv/2Lu), anit) ^ -012(1 - n/2uj), (11) 
Q!3i(t) = —a32{t — Tv/2uj), a[i2{t) — +031 (t — 7r/2ti;). 

From the first two lines we conclude firstly that an, Q12, 021 and 
Q22 have as functions of time a basic period of tt/ci; (not 2n/uj) 
and that a22(i) = aii(f±7r/2aj) and 0:21 (i) ~ "ai2(i ±7r/2(x;). 
The last line of dllb tells us that the averages of 031 and 032 over 
the period 27r/tt; vanish so that 031 and Q32 are simply oscillations 
around zero, and that they change their signs under time shifts by 
n/u). Our reasoning for aij applies analogously to rjij. 

We write down the result of these considerations in the form 



an 




022 


= a{t - 


- tt/2lo) 


ai2 - 


= ~m 


Q21 


= 7(i- 


- tt/2u;) 




= k{t) 


Q32 


= k{t - 


- 71/210) 


7711 




7/22 


= Vt{t- 


- %/2lo) 


V12 = 


= -m 


7/21 


= 5{t- 


■k/2uj) 




= A(t) 


V'i2 


= \{t- 


7r/2aj) 



(12) 



Here a, 7, fjt and 5 are in general periodic functions of time with 
the basic period ti/uj, but k and A are periodic functions with period 
2n/uj, which show sign changes under any time shift hy -k/lu and 
vanish under averaging over the period 2-k/uj. 
From (|9} and M2\ we conclude 

S = a{x ■ B)x + {y ■ B)y 
-j{y ■ B)x + 7^(5; ■ B)if 

-fjt [x ■ J)x - i)l {y ■ J)y (13) 
+S{y-J)x - ~S\x-l)y 

+ [k{x-B) + k\y-B) + ~X{x-l) + V{y-l)\ z. 

Here fi^, 7^, fjl and 5^ differ only by a phase shift of 7r/2 from a, 
7, fjt and 5, respectively, and and by a phase shift of tt from 
k and A. 

In addition to fields as defined above by averaging over x and 
y we consider also time-averaged mean fields defined by additional 
averaging over a time interval of length 2tt/uj (but we refer to them 
only if explicitly indicated). When speaking of time averaging in 
what follows we always refer to this interval. For time-averaged 
mean fields il3t turns into 



e = a[B - {z ■ B)z\ +-izx B- i-jtJ - 5z x J , 



(14) 



where a, 7, rjt and 5 are time averages of a, 7, fjt and 50 

In the special case of the Roberts flow, i.e. e = 0, the coeffi- 
cient a is independent of time and so coincides with a^, and this 
applies analogously to 7, fj, S, k and A. In addition in this case the 
inversion of i; in (|2l( is equivalent to a shift of the flow pattern, e.g., 
by tt/ V2 kn along y — x. Since such a shift does not change aver- 
ages, £ as given by J14l l must be even in z. Therefore we have then 
7 = 5 = 0. Nonzero 7 and 5 terms in J14t require a break of this 



^ In view of the signs of a and 7 we deviate here from representations as 
given, e.g., in Radler et al. (2002a) but follow CHT06. 



symmetry, that is, a preference of z over —z, and this may occur 
as a consequence of the aforementioned circular motion of the flow 
pattern. 

We override for a moment our restriction to non-negative val- 
ues of the frequency uj and admit also negative ones. For a; > 
the circular motion of the flow pattern defines, together with the z 
direction, a right-handed screw, and for < a left-handed one. 
We conclude from this fact that inversion of the sign of u has no 
other consequences than inversion of the signs of 7 and 5. 



Second— order approximation 

The task of determination of £ is now reduced to the determina- 
tion of the six functions a, 7, fjt, 5, k and A which occur in l ll2t 
and (Tsjl. As a first step in that direction we investigate £ within 
the second-order correlation approximation (SOCA). Later we will 
proceed to a corresponding fourth-order approximation. 

SOCA is defined by the neglect of the term with uxb — u x b 
on the right-hand side of equation ([8} for b, which turns so into 



(9t-77V^)6= (B- V)u-(u- V)B, V-6 = 0. 



(15) 



We may solve this equation with u as given by Q and {S} analyt- 
ically and calculate then £, see Appendix iBl When choosing the 
form ( 113b of the result we have 

a = ?io-RmX*^'(*) , 

Vt = iwofcn'-Rm [x*''W + x'"(i- V2a-)] , (16) 
7 = (5 = ft = A = 0. 
Here we have used the definition 

Rm = uo/rjkii, (17) 
and x'^' is given by 



where 



CCiujt , Lot — qr) e. ^ dr 



CC(a, fe) = cos [e(cos a — cos 6)] 
and 

q = Lo/rjk^ . 



(18) 



(19) 



(20) 



The parameter q gives, apart from a factor 27r, the ratio of the de- 
cay time of a magnetic structure with a length scale 2-n /kn, that is 
{2-k)'^ / Tjk^, and the wobble period 27r/aj of the flow pattern. In the 
case of small q the magnetic field follows the fluid motion imme- 
diately, but for large q it does so only with large delay. These two 
cases are sometimes labelled as "low conductivity limit" and "high 
conductivity limit", respectively. 

In agreement with the general findings summarized in M2\ . 
the function x'^' is periodic in time with a basic period h/lo. 
Whereas a and ct^ differ by a phase shift of 7r/2, fjt and fjl co- 
incide, x'^' satisfies |x'^^| < 1- It must be positive as long as 
e < 7r/4 but may otherwise take negative values, too. If e = 0, or 
e 7^ and g = (what corresponds to the low-conductivity limit), 



X 



(2) 



is independent of time and equal to unity. In Appendix iBlsome 
numerically determined values of x'^' are given. We note further 
that 



(2) 

X 

and 



(t) = 1 — {eq) sin ujt if eg <C 1 and q , 



(21) 
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x'^^(t) ^0 as q 



(22) 



For time-averaged mean fields we have again l ll4t . now with 

OL = uoRmxi'' , Vt = uok^^ Rn^xlP ; 7 = ^ = 0, (23) 

where Xo^' means the time average of x'^' over the period tv/lj. 
We point out that the time average of a function of cut, say f[uot), 
over an interval of the length n/u) is independent of lu. This is obvi- 
ous from {uj/tx) J^^^ f{ujt) dt = (I/tt) f{<p) dip. Hence x[P 
does not explicitly, but only in a indirect way via q, depend on u). 

(2) 

If e = 0, or e 7^ and g = 0, we have Xo — 1- For e = we 
fall back to the Roberts flow. Indeed, the result l l23t with Xo ^ ~ 1 
agrees with earlier results for this flow; see Appendix lAl In view of 
l l23t we note further 



'^'-l-|(eg)' if eg<landg<l 



X6 
and 



(2) n 

Xo ^ as 



(25) 



The last statement implies a/uoRm — > as g ^ oo. 

As we know from general considerations on SOCA (e.g., 
Krause & Radler 1980) the range of applicability of SOCA de- 
pends on q. For small q a sufficient condition for its validity reads 
Rm <^ 1. For large q such a condition is Rm/q <^ 1. 



Higher-order approximations 

Going now beyond SOCA we start again with Eq. ([8) for b and put 

b = + + + • • • (26) 

with b'"' being of the order n in u, and correspondingly 

= £ , = (uxb*"') , n> 1.(27) 

In that sense b and £ in Sect. l2.3l have to be interpreted as b'^' and 
f respectively. 

From ^ and l l26b we obtain l ll5t . now with 6*^^' instead of b, 
and further 

-V X (u X , (28) 
V-6'"+''=0, n>l. 

Using our result for and ^ we have calculated b'^'. It 
turns out that the average of u x b'^' vanishes, that is f '"^^ = 0. 
In the same way we may calculate b'''' and £^*\ However, these 
calculations are rather tedious. For the sake of simplicity we have 
ignored all contributions to S^'*'^ resulting from derivatives of B, 
that is, the terms with fjt, 5 and A in U3b . Some details of the cal- 
culations are explained in Appendix Icl 

Considering the results of all approximations up to the fourth 
order and referring again to l |13l l we have now 

5 = uoRm[x'^^Hi)~kRi^x'-*''\t)] , 

7 = ^uoRlx^^''\t), k = 0. (29) 

The functions x^* "'^ and x*^^ are given by 



p OO poo p oo 

2 CC(L0t,iut-q{r' + r" + t"')) 

Jo Jo Jo 



CS{LJt - qr' ,u;t - qir' + t")) 
exp(— (r + 2r -|- r ) dr dr dr , 



Ml) 



= 2 



OO /"OO /"OO 



^C{ujt,ujt — q{r' + r")) 



(30) 



SS{ujt — qr' ,Lot — q{T' + r" + r'")) 
exp( — (r' + 2r" + r'") dr' dr" dr'" , 

with CC as defined by ( 119b and analogously defined quantities CS, 
SC and SS, 



CS (a, 6) — cos [e(sin a — sin fe)] 
SC (a, 6) — sin [e(cos a — cos fe)] 
SS (a, b) — sin [e(sin a — sin b)] . 



(31) 



Note that CC and CS are symmetric but SC and SS antisymmetric 
in the two arguments. 

Like x''^' both x'*™' and x'*''' oscillate with a basic period 



(24) n/u. They satisfy |x'''°^i < 1 and jx'^'^'j < 1. Further x'''"' is 
positive as long as ]ej < 7r/4. In contrast to x'''"'*' however, the 
time average of x'*^' over a period n/ui is equal to zero. Whereas 



x''*"' is even, x'*^' is odd in lu. We have further 



and 



= 1- |(eg)^(l + 16sin''a;t), 
= ^(eq)^ sin ujt cos Lot , if eg ^ 1 and g^l, (32) 



,X 



(47) 



as g 



(33) 



For time-averaged mean fields again relation lll4b applies, 
now with 



(2) _ 1 p2 
AO 2 -"^mXo 



(4a) 



7 = 0, 



(34) 



where Xo 'i™^ average of x " . Like Xo also Xo " does 

not explicitly depend on ui. Unfortunately, values for rjt and S are 
not available. We have 



X^*^' = , if eg < 1 and g < 1 , (35) 



(4q) (47) 
AO 1 AO 



and 

as g — > oo . (36) 
With l[34) we find then 

a = UoRn. [l - |(eg)' - - |(eg)2)] 

if eg < 1 and g < 1 . (37) 

Results of higher approximations are very desirable but re- 
quire heavy efforts. We suspect that in the approximation of sixth 
order in u the time averages of 7 and 5, and so the coefficients 7 
and S in JI4b no longer vanishes. This presumption is supported by 
numerical results (see below). 



2.4 Mean electromotive force in case (ii) 

Modifying the considerations on case (i) properly we may conclude 
that relation ( 114b . again considered for time-averaged fields, applies 
for the fluid flow of type (ii) with 7 = 5 = 0. By contrast to case 
(i) the correlation between velocity components at different times 
vanishes if the time difference becomes very large. 

Modifying also the SOCA calculations described above and 
in Appendix |B] correspondingly we find again ( 123b . but with Xo^' 
being the time average of 

/•oo 

X^^\t)= / cos{e[04t/rc)-04t/rc-gr)]}e'"dr, (38) 
Jo 
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Figure 1. Dependence of Xq on g for two different e, calculated numer- 
ically on the basis of equation j38t . The dashed line shows that Xq"' with 
e = vr behaves like q^^ for large q. 



where q is now defined by 

q={r,rikl)-\ (39) 

We have here again x'^' = 1 for q = 0, and x'^' vanishes for q 
oo. Like x also Xo depends on e and g but no longer explicitly 
on Tc. We have calculated Xo^' on the basis of equation l|38) under 
the assumption that (j)^ is always constant over time intervals of a 
given length. Fig.[T]shows dependencies on e and q. 

3 TEST-FIELD METHOD 

We will determine numerically the elements of the tensors a^j and 
r]ij introduced with l|9), but with 1 < i, .7 < 2 on ly, employing the 
test-field method o f ISchrinner et"ai] ( l2005ll2007t) . 

We will calculate £ — u x b from numerical solutions b of 
([8}, with B replaced by one out of four test fields B^', 

B^'^ = B {cos kz, 0, 0) , B^" = B {0,co8kz,0) , 

B^' = B {sin kz, 0,0), B^' = B (0, sin fcz, 0) , (40) 



where B and k are a constants. Repeating this for all test fields, 
denoting the £ that belongs to a given B'"' by f', and using {9} 
we find 



(41) 



(42) 



^Ti^) ^ B {'^ip cos kz ~ r^Jpfc sin kzj 
^Ti^) ~ ^ (o^ip sin fcz: + ry^^pfc cos kz) 
for 1 < j, J < 2, where 

t (-1112 VU 

\^-r/22 7?21 

From this we conclude 

Qij = B^^ ^^^'^{z) COS kz + £i^{z) sin kz'j , 

vlj = -{kBy^ \£l''{z)sinkz-£l\z)coskz] , (43) 

again for 1 < j,i < 2. 

We point out that, although the £'"' depend on z, the aij 
and rjij have to be independent of z. We further note that re- 
lation on which these considerations are based, can only be 
justified under the assumption that all higher than first-order spa- 
tial derivatives of B are negligible. The derivatives of order n of 



our test fields B are proportional to fc". For this reason the re- 
sults l |43l l apply in a strict sense only in the limit fc ^ (cf. 
[Brandenburg, Radler & Schrinner 2008). 

Let us focus here on case (i). After having calculated the Uij 
and rjij in the way indicated above we may determine the a, 7, fjt 
and 5 according to l ll2t . that is, 

a{t) = i [qii {t) + a22 {t + 7r/2a;)] , 
7(t) = -\ [qi2 {t) - 021 {t + 7r/2a;)] , 

flt{t) = \[llll{t)+n22{t + 'K/2L0)] , (44) 

hi) = -^[rii2{t)-ri2i{t + Ti/2u)] . 

We are, however, mainly interested in the time-independent coeffi- 
cients a, 7, 77 and 5 that are relevant for time-averaged mean fields 
as addressed in l ll4t . They are just time averages of the ct, 7, f]t and 
5, that is 

a = |({qii) + (022)) , 7 = —^((012) — (021)) 

T), = i({r?ii) + {7;22)), 5 = -\{{-ni2) - imi)) , (45) 

where {■ ■ ■) means averaging over a time interval of length -n /uj.ln 
case (ii) the relations l |45t apply with (■ ■ •) interpreted as averaging 
over a sufficiently long time. 



4 A GENERALIZATION 

So far we have assumed that the mean electromotive force £ in 
a given point is completely determined by B and its first spatial 
derivatives in this point. If we relax this assumptio n we may pro- 
ceed as in 'Brandenburg, Radler & SchrinneJ boosi) . In that sense 
we may replace l|9ll, applied to time-averaged mean fields, by 

£,{z)= j (C)B, {z-O-f,,, (C) J, [z-O] dC (46) 

with kernels &ij and fjij. When using a Fourier transformation 
Q{z) ~ J Q exp(ifc2)d2;, this turns into 



£,{k) ^ &,,{k)Bj{k) - fj,j{k)J,{k) , 
where 



(47) 



Qij (fc) =y (0 ^'^^ ^( , rjij (k) =y ''^^ (0 COS k( d( . (48) 

In this understanding the relations (l4U-ll43t apply with aij and rjij 
being replaced by Oij and fjij , which have a well-defined meaning 
for all k (not only in the limit k —> 00). 



5 RESULTS 

5.1 Units and dimensionless parameters 

It is appropriate to give a and 7 as well as a and 7 in units of uo, 
and fjt, S, rjt and 5 in units of uo /kn- The remaining dimensionless 
parts of these coefficients are then, apart from the time dependen- 
cies of a, 7, fjt and S, functions of the dimensionless parameters 
-Rm, e and q introduced through J17t . lO, and either l l20t or l |39t . 
Instead of q we may also use the dimensionless quantity Cj defined 
by 



Cj = q/Ru 



(49) 



© 2008 RAS, MNRAS OOO.mflOl 



6 K. -H. Rdcller and A. Brandenburg 




jg 2Q q = 0.1. The dotted curve coiTesponds to 137) , which has been derived for 

A^ eg < 1 and g < 1 only. 



Figure 2. Time dependence of a^j for the parameters used in Fig. I of 
CHT06 which, in our normalization, are i?m = 78, e = 3/4, and Cj = 
^2/3. Here, At = t — to> where to = 300/w is the final time shown in 
Fig. 1 of CHT06. The dotted lines refer to an and a2l. respectively, the 
dashed lines to a\2 and 022, and the dash-dotted lines to (an + a22)/2 
and (021 — oi2)/2. The straight solid lines give the time averages of the 
latter quantities, that is, a and 7. 



In case (i) we have so a; = u/uoku, which is the ratio of the 
turnover time (uokji/^-K)^^ to the wobble period 27r/aj. In case 
(ii) apphes a) = (rcUoku)^^ , and this is, apart from a factor 2tt, 
the ratio of that turnover time to the time Tc introduced with the 
random flow. 



5.2 Case (i) 

Comparison with CHT06 

We show first that our method reproduces results by CHT06. We 
suppose that our Rm is related to the magnetic Reynolds num- 
ber, say R^^^ , used but not explicitly defined there, by Rm = 
While in CHT06 dependencies of the results on Rm 
and e are considered, no values of g or d> are given. We suppose that 
the calculations have actually been carried out with uj = -^2/3. Fi- 
nally we suppose that the unit of a and 7 used by CHT06 isuj/ku- 
With a view to Fig. 1 of CIIT06 we have carried out calcula- 
tions with Rm = ^/3/2 X 64 « 78, e = 3/4 and li = y/VS. Our 
results for the Uij obtained with these parameters and given in this 
particular case in units ofui/ku are presented in our Fig.|2] We see 
in particular that an and 022 vary between —6 and —1 with a pe- 
riod n/u!. As far as qh is concerned this agrees with the result for 
{u X b)x shown in Fig. 1 of CHT06. Also the initial evolution of 
an, which is not shown here, agrees with this figure. Furthermore, 
in our Fig.[2]the phase shift by n/2 between an and 022 discussed 
in Sect. l2.2l is clearly visible. Our results for Qfi2 and 021 lead to a 
value of 7, which agrees in modulus but differs in sign from that of 
CHT06. (To obtain their sign we need to replace lj by —uj.) With 
the above values of Rm and lu but e = 1 we find again a sign of 7 
opposite to that of CHT06. 



Time-averaged mean fields 

Switching now to time-averaged mean fields we start with Fig. [3] 
which shows results for a at Rm = 0.1 in dependence on eq. They 
were found with the help of numerical integrations of the test-field 
version of ((Sj in its complete form or after reducing it to SOCA. 
It turned out that SOCA is sufficient for their calculation. Some of 
these results were also confirmed by evaluating l l23b with l llSt or 
i24l . As long as eq is small, a depends in agreement with i23l and 
l l24b only via this product on e and q. For larger eq it depends, how- 
ever, in a more complex way on e and q. Furthermore, a remains 
finite if e = 1 and q grows, and it tends to zero if g — 0.1 and e 
grows. Since Rm is small the validity of SOCA is plausible in the 
case g = 0.1. It is however remarkable in that with e = 1, in which 
g may grow up to 10. 

Next, we consider the dependence of a and 7 on Rm, in Fig.|4] 
shown for e = 1 and Cj = 1 (i.e. g = Rm)- For small Rm we expect 
that SOCA applies and so a/uQ is linear in Rm but 7 vanishes. In- 
deed a/uit shows this linearity up to Rm ~ 1. In agreement with 
the results of CHT06 7 is negative and its modulus remains small 
for Rm < 1. Remarkably the values of a/uo calculated from ( 1231 ) 
and J24b (dotted line), or ([37} (dashed line), which have been de- 
rived for g ^ 1 and eg <C 1, deviate for Rm > 1 drastically from 
both the numerically obtained SOCA results (dash-dotted line) and 
those obtained without any approximation of that kind (solid line). 
The proportionality of 7/110 with Rm confirms the presumption 
made at the end of Sect. l2.3l that nonzero values of 7 occur only in 
sixth-order and higher approximations with respect to ua- 

Simple arguments (as given in Sect. |6] below) suggest that a 
is never negative. However, CHT06 found that not only the mod- 
uli but also the signs of both a and 7 depend for each given Rm 
sensitively on e. In our Fig. [5] which applies for Rm = 100 and 
Co = •\/2/3, both a and 7 vary strongly with e, too. The repre- 
sented results confirm, apart from the sign of 7, the corresponding 
ones in Fig. 2 of CHT06. Both a and 7 change their signs with e. 
As Fig.|6]shows, in the situation with the same Rm and lj — 1 only 
7 changes its sign, which indicates a considerable effect of chang- 
ing dj. In both of the cases considered in Fig.|5]and Fig.|6] a and 7 
diminish for small as well as large values of e. 

In Fig. Hand Fig. [8] we see that a and 7 depend, at least for 
Rm ~ 100 and e = 1, also sensitively on the parameter ti, or g, 
that is, on the frequency with which the velocity pattern wobbles. 
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Figure 4. Dependence of a/no and —■y/uQ on Rui for e = 1 and a; = 1 
(i.e. q = Rm)- Solid lines show results obtained without any approxima- 
tion. The dash-dotted line gives a/uQUS obtained numerically using SOCA. 
The dotted and the dashed line give results calculated with i23\ and j24t . or 
)37t , respectively. 
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Figure 5. Dependence of a/uo and —■y/uo on e for Rm = 100, and 
ill = \/2/3 (the value considered by CHT06). 



Figure 7. Dependence of a/ug and —■j/uq on a; for i?m = 100 and 
e = 1. 
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Figure 8. Logarithmic representation of the dependence of +a/uo (thick 
solid lines), —a/uQ (thin solid lines), —"f/uQ (thick dashed lines), and 
+7/jiO (thin dashed lines) on ui for Rm = 100 and e = 1. For large 
values of ui we have a/uo ~ 0.065 and —-y/uo f» 0.8/ui. 




Figure 6. Dependence of a/uo and —y/uo on e for Rm = 100, and 
li = 1. 



There are, however, simple asymptotic behaviors for small and for 
large ui, clearly visible for < 0.1 and ui > 3. Similar results have 
been found for Rm = 10 and e = 1. In this case, however, a stays 
positive for all values of ui, and only one sign reversal of 7 occurs. 

We see from CHT06 that there is a rich dependence of a and 7 
on Rm for values of ti and e of order unity. In Fig.|9]we show results 
for an example with Cj = 0.5. Reversals of a are then possible for 
rather small values of Rm of the order of 10. However, as Fig. 1 101 
shows, such behaviour disappears for Cj = 10, in which case a 
stays always positive and 7 always negative. In fact, there is an 
asymptotic scaling a/uo ~ Rm^^^ as Rm 00, and 7 approaches 
a constant finite value as Rm — * 00. 

In a few cases rjt and S have been determined in addition to 
a and 7. Results on the dependence of these quantities with e — 1 
and Lj — 0.7 on Rm are shown in Fig.[TT] They have however been 
calculated with k = fen, not fc ^ 0, and are therefore at most 
approximations of the mentioned quantities. 

A correct interpretation of these results requires a look on the 
explanations of Sect.|4]on the non-local connection between £., B 
and J as defined by i46b . In that sense the a, 7, rit and 5 in Fig.llll 
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line), —a (thin solid line), —7 (thick dashed line), and +7 (thin dashed 
line) on /?m for tj = 0.5 for e = 1. 
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Figure 10. Dependence of a and 7 on Rm for a) = 10 for e = 1. Note 

— 1/2 

the asymptotic scaling a/uo ~ (dash-dotted line) and that 7 ap- 

proaches a constant finite value as Rm — ► 00. 

may be understood as values of the functions a{k),-Y{k),%{k) and 
S{k) at fc = kn- 

In the following, when writing a{k) or rit{k), for example, 
we always mean a{k) or fjt{k). In an earlier investigation with 
the Roberts flow u nder SOCA and with isotropic turbulenc e inde- 
pendent of SOCA terandenburg. Radler & SchrinnedllOOSi) It was 
found that a{k) and rit{k) vary with k i n a Lore ntzian fashion like 
( 1+fc^ ) ' \ However, for the GP flow lCourv oisier ( 2008) found 
that a{k) at small k is extremely sensitive to the value of Rm- 

Fig. 1121 shows that a{k) and "/{k) for Rm = 30, Cj = 0.5, 
and e = 1, approach the values given in Fig.|9]as fc ^ 0. However, 
the magnitudes of rjt (fc) and S{k) become rather large as fc — > 0. It 
turns out that a is positive for fc/fcn > 0.4 and 7 becomes smaller 
with increasing fc. Remarkably, J]t(fc) is negative for fc/fcn < 1, 
suggesting that magnetic field generation might be possible via a 
negative magnetic diffusion instability. 

In order to check this possibility we have calculated the linear 
growth rates 

A±(fc) = -[r7 + r?t(fc)]fc'±«(fc)fc. (50) 
Fig.|13|shows that A± is almost entirely given by ±Q(fc)fc. A neg- 



Figure 11. Dependence of a and 7 (dashed lines) as well as rjt and S (solid 
lines) on Rm, for lu = 0.7, e = 1, and k = kn. 
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Figure 12. Functions a(fc), 7(fc), »;t(fc), and (5(fc) for iJm = 30, ui = 0.5, 
and e = 1. 



ative diffusivity instability does not occur It is important to realize 
that most of the small wavenumber modes, especially those with 
negative values of rjt, would never be realized. This is because in a 
system of given size, only the corresponding harmonics will have 
a chance to be excited, and of those only the ones with the largest 
growth rates will dominate. We should point out that the detailed 
variations of A-|- shown in Fig.[T3]may not be accurate. In fact, this 
figure shows a maximum at fc/fcn ~ 0.75, but direct simulations 
suggest that the fastest growth occurs for fc/fcn ~ 0.5 with a grow 
rate of A ~ 0.23urmsfcf ■ Nevertheless, this value is still compatible 
with Fig. [T3] 
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Figure 13. Dependence of \± on k for Rn 
The dotted lines give ±ofc for comparison. 
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Figure 14. Dependence of a and 7 on li for the flow with random time 
dependence (case ii), with i?ni = 100 (solid lines), iJni = 10 (dashed 
lines), = 1 (dash-dotted lines), and R^n =0.1 (dotted Unes), and 
with e = TT. The error bars are similar in all cases, but are only shown for 
i?m = 100. 



5.3 Case (ii) 

In case (ii) we have calculated a and 7 under the assumptions on 
4>x and 4>y introduced in Sect. 12.41 Fig. [14] shows results for Rm 
ranging from 0.1 to 100 and e = tt as functions of Co. In the limit 
of small Co the flow can be considered as stationary, that is, as a 
Roberts flow. Indeed in this limit the va lues of a agree well with 
those obtained for the Roberts flow (e.g.. lRadler et alj ( l2002ah . see 
also Appendix IaI. For large values of Co the values of a vanish for 
all Rm. For not too small Co and _Rm there is no longer a noticeable 
variation of a with _Rni, and a reaches a maximum at (D ~ 0.3. 
In the range 0.3 < Co < 1 continuous flow renewal removes the 
tendency for a to diminish with growing R^. For Co < 0.2 the 
value of a remains strongly dependent on Rm and can still change 
sign. We have also calculated a with e — 1 and Rm ~ 100 as 
a function of uj and found a qualitatively similar behavior as for 
e = TT. In this case it remains positive and is up to 50% smaller 
than for e = TT when Co < 1 and somewhat larger when Co > 1. In 
all cases we found, as expected, 7 = within error margins. 



6 DISCUSSION 

Our results for the flow of type (i) confirm the finding of CHT06 
that both the a and 7 coefficients depend sensitively on Rm and 
also on e, and that even the signs of these coefficients may vary with 
these parameters. We have to add that a and 7 depend also on Cj, or 
q — Cj Rm, that is, on parameters connected with the frequency of 
the wobbling motion, which CHT06 fixed in a special way without 
commenting on it, and that they show similar variations with these 
parameters. We found however rather regular behaviors of a and 7 
for small and for large values of e and Co. 

It is sometimes considered as a rule that the sign of a is op- 
posite to that the mean helicity of the fluid flow and its modulus 
is proportional to that of the mean helicity. There is however no 
general reason for that kind of relation between a and the kinetic 
helicity. We see only two limiting cases which allow simple state- 
ments on the sign of a. 

Firstly, in the framework of SOCA applied to homogeneous 
turbulence and comparable flows it turns out that the sign of a in 
the limit q — > is always opposite to th at of ■!/) ■ (V x ■0) , wher e 
u = V X V ■ = 0; see, e.g.. iKrause & Radlerl ( IT980I) : 
iRadler & Brjuidenburj Ewf). For both types of flows, (i) and (ii), 
we have if) = {ip, —tp, %p), where ip = —kii(d^/dx + dip/dy) 
and therefore 1/) • (V x 1/)) = —2u^/ku_. This implies that a is 
positive. Indeed, only positive a have been observed for small g, 
even beyond SOCA. 

Secondly it was found in SOCA under the same 
conditions, but in the limit q — > 00, that the sign of 
a for a flow with finite correlation time is opposite to 
that of r u(x,t) ■ X u(x,t ~ t) dr; see, e.g., again 
Krause & Radleii ( 119801) . A relation of that kind between a and 
this integral ca n indeed be formally derive d from the general 
relation (29) of iRadler & Rheinhardtl ( l2007h and applied to our 
specific situation. In case (i) the correlation time is however 
infinite and this integral does not converge. Although we know that 
w ■ (V X u) — — 2itofcH we do not see how reliable conclusions 
could be drawn concerning the sign of a in the limit q ~* 00. In 
case (ii) the integral is positive, and indeed only positive a have 
been observed. 

Beyond the low and high conductivity limits, that is, for not 
too small or not too large values of q, even SOCA offers no simple 
general statements on the sign of a. In general a may take both 
positive and negative values. 

For studying a and rjt with very simple flows it seems appro- 
priate to consider flows of type (ii) rather than of type (i). In case 
(i) the results are influenced by the aforementioned circular motion 
of the flow pattern. As long as only a and rjt should be discussed 
there is hardly a reason to introduce such a motion. We see no nat- 
ural interpretation of it and so no interpretation of the so caused 7 
and 5 effects. 

Recently ^T ilgned ( 1200 8 l) pointed out that a time-dependent 
flow of a conducting fluid can act as a dynamo even when steady 
flows which coincide with it at any particular time cannot. He 
demonstrated this with a Roberts flow modified by a drift of its 
pattern so that the velocity u satisfies relations like ^ and ([3) with 
ffix = —knv^t and tpy — 0, where is constant the drift ve- 
locity. Even if the intensity of the flow is too weak so that in the 
case Wd = no growing solutions of the induction equation with 
a given period in the z direction exist, such solutions may occur 
in an interval of some finite v^. Although this flow considered by 
Tilgner is in a sense simpler than the flows in our paper, it shows 
no longer the symmetries with respect to the z axis which we have 
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utilized. As a consequence the relation between the mean electro- 
motive force £ and the mean magnetic field B is more complex. In 
particular l |13l l and il4t no longer apply. Nevertheless the question 
arises whether the effect of the time-dependence of flows observed 
by Tilgner occurs also in the examples investigated here. In case (i) 
the parameter uj could play the role of . The fact that the magni- 
tude of a is larger for some finite Cj than for d; = 0, which can be 
seen in Figs|7]and[8] points in this direction. 

One of the original motivations for looking at the GP flow 
was the fact that it is time-dependent and in that sense closer to 
turbulent flows than time-independent flows. However, as we have 
shown here, the dynamo properties of the GP flow cannot be com- 
pared in a meaningful way with analytic theories or with simula- 
tions that apply to isotropic turbulence. Nevertheless, as shown in 
this paper, an analytic theory for the a effect and other turbulent 
transport coefficients can be derived that matches numerical results 
in limiting cases. 

Astrophysical flows can often neither be described by 
isotropic turbulence nor by wobbling two-dimensional flow pat- 
terns, but they are likely to contain aspects of both extremes. 
However, the present work highlights another aspect that may be 
of more general significance and concerns the turbulent transport 
properties in the presence of high-frequency time variability. This 
is not just a peripheral aspect of turbulence, but it is an additional 
property whose effects need to be understood more thoroughly. The 
situation is reminiscent of the modifications of mix ing length the - 
ory in the presence of stellar pulsations (see, e.g., lGouglj|l977h . 
In dynamo theory the issue of high-frequency time variability has 
only recently been addressed. One example concerns the nonlinear 
a effect where its time dependence has a striking effect on the be- 
haviour of the mean field. In that example the temporal behaviour 
of the forcing function (delta-correlated or steady) determines the 
nonlinear asymptotic scaling behavio r of the quenc hing func tion 
a(B) at low R^. The early results of lMoffattI ( Il972 ) and RudigeJ 
ill 974) sug gested a |a| ~ \B\^ '^ behavior, but in more recent years 
iField et aL (1999) and Rogachevskii & Kleeorid ( I200Q) found in- 
stead a |a| ~ behavior, whic h seemed i n confl ict with the 
earlier results. However, the work of ISur et al] j2007 h now shows 
that this is not just an artifact related to different approximations, 
for example, but it depends on whether or not the flow is time- 
dependent. They found that the |a| ^ IS]^"^ behavior is repro- 
duced if the flow is steady, while the \a\ ~ |-B|^^ behavior is 
obtained in the time-dependent case using a forcing function that is 
(5-correlated in time. Again, it is not clear which types of flows are 
more astrophysically relevant, but it is now clear that the detailed 
time-dependence of the turbulent flows can affect its transport prop- 
erties in rather unexpected ways. 
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APPENDIX A: ROBERTS FLOW 

In the special case e = the Galloway-Proctor flow, defined by ^ 
and l[3), turns into the Roberts flow. Our SOCA results for this spe- 
cial case agree with results for the Roberts flow reported in Bran- 
denburg et al. (2008) and in Radler et al. (2002a), referred to as 
BRS08 and R02a, respectively. 

In BRS08 instead of our coordinate system (x, y, z) another 
one, say {x' , y' , z), is used, which is obtained by a 45° rotation of 
our system about the z axis, that is, 

X ^ -^{x - y') , y = ^{x' + y'). (Al) 

In the case e = 0, to which we restrict ourselves here, {5) turns 
under this transformation into 

t/j = ^ cos {kjix'/V2) cos {kHx'/V2) . (A2) 

Together with l|2j we find so, referring to the system {x' , y' ,z), 

I -cos{kHx'/V2) sm{kuy'/V2) \ 
u — V2uo I +sin cos (kny' . (A3) 

V-\/2cos {knx'/V2) cos {kHy'/V2) j 

Comparing this first with (BRS08 25) and ignoring the oppo- 
site sign of Uz we find 
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BRS 



ko = kf/V2 = kH/V2. 



(A4) 



The only consequence of inverting the sign of is a sign change 
of a. Taking then the SOCA results (BRS08 30) for a and r]t, with 
itQ in place of ito and completed by (BRS08 29), considering 
l lA4t and the remark on the sign of a we can easily reproduce our 
results l l23t . 

This applies analogously to R02a if, in addition to the trans- 
formation l lAlb . x' is replaced by x' — 7r/\/2fcH and y' by y' + 
7r/-v/2fcH. Comparing the corresponding modification of l lASt with 
(R02a 15) we find 

u± = {2V2/tt)uo , U|| = (8/7r^)uo , a = \/27r/fcH . (A5) 

When using (R02a 19) we obtain our result ^ for a. With (R02a 
38) and rjt = l3± + we may also reproduce our result J23b for 
Vt- 

Going beyond SOCA we note that according to (BRS08 25), 
or also according to (R02a 20), 

Q = UoRm(l}{2Rm) (A6) 

with a function (j> satisfying <j!>(0) = 1 and vanishing like R^^^^ 
with growing _Rni- It has been calculated numerically and is plotted, 
e.g., in R02a. 



APPENDIX B: SECOND-ORDER CALCULATIONS 

For the calculation of the Q^j and r]ij with 1 < i, j < 2 under 
SOCA we start with M5\ . Introducing there b = Re(f) exp ifcz) 
and B = Re(i? exp ikz) we obtain 

[dt - ??(V^ - k^)] b={B ■ V)u - -iku^B . (Bl) 
For our purposes it is useful to represent u in the form 

Ux = -ito[cos(fcHy) ss(t) +sin(fcHy)cs(t)] , 
Uy = Mo [ cos(fcH3::) sc(t) + sin(fcHa;) cc(t)] , 
Uz = — Mo [ cos(A:hx) cc(t) — sin(fcHa;) sc(t) (B2) 
+ cos(fcHy) cs(t) - sin(fcHy) ss(t)] , 

where 

ss(f) = sin(e sin Ljt) , cs(t) = cos(esinLjt) , etc. (B3) 

Then the right-hand side of JBlt . say i?, takes then the form 

R = H cos knx + H sm kh^:: 

+R. " cos kuy + sin kuy (B4) 

with .R , R , ■ ■ ■ depending on time. 

Clearly, JBU poses an initial value problem. As initial time 
to we take to — oo. Then the solution b of JBlt is completely 
determined by its right-hand side, R, and has again the form of R 
as given by l lB4t . Since — —k^b we have 



bit) 



R{t - t') exp[-r7(fci[ + At' . 



(B5) 



After determining the R , R , 

" cx "ax 

defined b , o , ■ ■ 



we find for the analogously 



blj^ = iMofc-Ba;r(cc) , b^x ~ ~iuoki3xT{sc) , 

bl^ = —uo{kiiBy — ikBx)T{cs) , 

b^x = uo{kjiBy — ikBx)T{ss) , 

b^y = uo{kjiBy +ikBx)T{cc) , 



Table Bl. Some values of x'^^ defined in equation (Ts). 



TV / O 


e = 0.1 


e = 0.1 


e = 1 


e = 1 


e = 10 


e = 10 


9 = 1 


<2 = 10 


q = l 


(J = 10 


q = l 


g = 10 





0.9970 


0.9926 


0.7558 


0.4204 


0.3017 


0.2073 


1 


0.9982 


0.9936 


0.8521 


0.4882 


0.5246 


0.2667 


2 


0.9990 


0.9954 


0.9123 


0.6165 - 


0.1882 


-0.1816 


3 


0.9989 


0.9971 


0.8963 


0.7341 


0.1567 


0.1912 


4 


0.9980 


0.9975 


0.8115 


0.7658 — 


0.1558 


— U.2440 


5 


0.9968 


0.9965 


0.7121 


0.6889 


0.1196 


0.1881 


6 


0.9960 


0.9947 


0.6586 


0.5550 - 


0.0363 


-0.1637 


i 


0.9961 


0.9931 


0.6775 


U.4400 


0.1317 


0.2189 


8 


0.9970 


0.9926 


0.7558 


0.4204 


0.3017 


0.2073 


by = 


—Ito 


{kuBy + ikBx)r{: 


5C) , 




(B6) 


b"^ - 

Uy 


iuokByFics) , 


^sy - 

Uy 


—iuokByT( 


>s), 




bT = 


uokiiBxT{sc) , 


tax 

bz = 


uokuBxri 


cc) , 




b? = 


uokuByl'iss) , 




uoknByFl' 


cs) , 




where 














r(/) = 


/ fit 

Jo 


— t') exp[ 


-Vikii - 


f k^)t'] dt' 




(B7) 



with any function / = f{t). Of course, F depends in general on 
time. In view of JB6b we note that, since B = Re(i3 exp ifcz), 
we have also J 

Jx- 

Calculating then £ we find 



Re( J exp ikz) and therefore ikBx = Jy and 



ikBy — 



£' X Uybz 



Uzby — UoRmaxxBx — UoRmbxxJx 



Uzbx — Uxbz 



£z 



Uybx 

■I , bxx 



UoRmClyyBy 

0. 



(B8) 



where axx, 
time, which are defined by 

axx = r]k^ [ccr(cc) + scF(sc)] 
ayy = r]ku [csr(cs) + ssr(ss)] , 



(B9) 



The combination of trigonometric functions on the right-hand side 
of the relation for Uxx can easily be expressed by the function CC 
defined by l |19l l. The same applies to ayy and the function CS de- 
fined in l l31t . 

The result given by l |B8t and l |B9l ) is valid for arbitrary k. 
This applies of course also if it is written in the alternative form 
with CC and CS. In that sense it is of some interest in view of 
the nonlocal connection between £ and B studied in the paper by 
[Brandenburg. Radler & Schrinneil ( l2008i") . In the main part of the 
present paper we consider however the limit A: ^ only. In this 
limit JB7b applies with fc = 0. Then ( IB8t and ( IB9t agree just with 
^ and OU. 



APPENDIX C: HIGHER-ORDER CALCULATIONS 

For the sake of simplicity we assume now, beyond SOCA, that B 
is a uniform field, that is, has no spatial derivatives. Then u x b is 
independent of space coordinates and ([28} turns into 

(n + l) _ f^{n) _ ^ V)b'"' 



(a 



rjV'')b'' 



{b^"> ■V)u■ 



n>l. (CI) 
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We may apply some modification of the procedure used in Ap- 
pendix |B]for solving the equation ( IB It for b to the equations l ICU 
for 6(2) and b*^). 

The right-hand side of the equation for 6*^^^ say R^'^\ is a 
linear combination of products ip{x) ip{y), where ip{x) stands for 
cos kux or sin knx, and ip{y) for cos kuy or sin kny. Clearly b'-^) 
has the same form as R^^^ and satisfies V^b'^-* = —2k^b^^\ 
Therefore l |B5l l applies after replacing b and R by b'^' and R^^\ 
respectively, fc^ by 2k'^, and putting fc = 0. As a consequence of 
the described structure of 6*^^' we have u x b'^' = 0. 

The right-hand side of the equation for b'"^', which we call 
fl'^', is a linear combination of products ipi{x) (p2{x) ifi{y) or 
ifi{x) ifii{y) ip2{y), where the indices 1 and 2 may refer to the 
same function or to different functions, e.g., ifi{x) — if2{x) = 
cosknx, or — cosknx and fiix) — sinkux. In 

the first case we utilize cos^ knx — |(1 + cos2A;Ha;) and 
split, e.g., cos^ knx sin kny into the two parts i sin kny and 
icos2A;HK sinfcny. In this way we may split l?'^-* into two 
parts, it*^"' and i?*^**', where /i*^"' contains only contributions 
ifiiix) ip2{x)ip{y) and v?(a;) ipi(y)ip2{y) with three different fac- 
tors, and contributions of the types '^{x) ip{2y) and if{2x) if(y), 
and ii'^''' only contributions of the types ip{x) and ip{y). There are 
two corresponding parts of b'^^^ , that is b'^"' and b'^*"' , which sat- 
isfy V^b*^"' = -5fc^b(="') and V^bf^*"' = -fc^ b''^") and e qua- 
tions of type of l lBSt . The structure of b'^"' impUes u x b*^''' = 
Only b'^*"', for which l |B5l l applies with b and A replaced by b'^*"' 
and R^''^''\ respectively, and A; = 0, contributes to u x 

Detailed calculations along these lines deliver us 



S'^^ = ""yfe?* — u^b^y^ = uo-Rm (a^2,_B^ +axyBy] 



Sy' = li^&^c — Uxbz = uoRfn[ayxBx + UyyBx'^ (C2) 

f i"** = - Uybl^^ = . 

with 

a^x = — (?7fc|)^{sc[r(ss,ss,cc) + r(cs,cs,sc)] 

+cc(r(ss, ss, cc) + r(cs, cs, cc)] } 

Clxy = — (?7fcH)^|sc[r(sS, cc, cs) — r(cs, CC, ss)] 

— cc[r(ss, sc, cs) — r(cs, sc, ss)] } 
ayx = +(77fc|)^|cs[r(sc,ss, cc) — r(cc,ss,sc)] (C3) 

— cc[r(sc, cs, cc) — r(cc, cs, sc)] } 
'^yy ~ — (?yfcH)'' |cs[r(sc, sc, cs) + r(cc, cc, cs)] 

+ss[r(sc, sc, ss) + r(cc, cc, ss)] } 

where 

nf,9,h) = 

f{t - t')g{t - t' - t")h{t -t' - t" - t'") (C4) 



OO POO POC 



exp[-rjkti{t' + 2t" + t'")] dt' dt"dt"' . 

The combinations of trigonometric functions in iC3l can be ex- 
pressed by the CC, CS, SC and SS defined in ([Till and ([IB- In this 
way we arrive at the results ([29} and i30\ . 



